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LAYER 


Donald 


SYNOPSIS 


new approach taken the problem incompressible, turbulent bound- 
ary layers. From physical considerations involving the diffusion turbulent 
stresses, the boundary layer divided into three regions distinguished the 
relative proximity the wall. Within distance about ten per cent the 
boundary layer thickness from the wall, the local wall conditions have im- 
portant effect; whereas the outer seventy-five per cent governed entirely 
the spatial history the boundary layer. The inner region includes the lami- 
nar sublayer, the buffer zone, and that part the turbulent region character- 
ized the logarithmic velocity profile, The outer region the turbulent 
boundary layer similar that wake jet, being characterized 
sharp but irregular boundary between turbulent and non-turbulent regions. 
Here the velocity profile may usually approximated three-halves power 
deficiency relation. The region between the inner and outer parts blending 
region more complicated than either the other two and not amenable sim- 
ple analysis. This new physical analysis, combined with empirical correla- 


tions, leads new method for calculating two-dimensional turbulent bound- 
ary layers that usually only involves algebraic equations and more versatile 
than any the previously published methods. 


INTRODUCTION 


The basic differential equations for turbulent flows were derived Osborne 
1895 assuming that all instantaneous velocities could ex- 
pressed the sums time-average velocities and instantaneous fluctuations 
about the average. Reynolds inserted his expressions for instantaneous veloc- 
ities and pressures the Navier-Stokes differential equations and took time 
average. The resultant equations include additional terms, known the Rey- 
nolds stresses, that involve the correlations the fluctuation velocities and 
which represent the turbulent transport momentum. When the Prandtl 
boundary-layer are applied these equations some the re- 
maining terms still involve the Reynolds stresses. because these terms 
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that mathematical solutions turbulent-boundary-layer problems have not 
been possible, they have been for the corresponding laminar flows. 

the past thirty years, various approaches have been taken turbulent- 
boundary-layer problems, all them without real success. The present sta- 
tus the problem that the theoreticians are awaiting the results hot-wire 
measurements turbulent shear flows with the hope that these experiments 
will reveal the basic mechanism turbulent stresses and thus enable the 
building rigorous theory. believed most workers the field that 
will many years before such theory completely developed and still 
more years before can yield results practical value the engineer. 

the meanwhile, engineering problems require answers and cannot wait 
for the ultimate theoretical solution, desirable that improved semi- 
empirical methods calculating turbulent boundary layers developed. The 
present paper constitutes report recent analytical study turbulent 
boundary layers aimed satisfying this this study, which was par- 
tially supported the Office Naval Research, the wealth experimental 
data and theoretical ideas the literature were reexamined and reevaluated 
terms some the newer information and ideas concerning turbulent 
shear flows. One result was the development new method for predicting 
the behavior two-dimensional boundary layers that applicable whenever 
the initial spatial variation the boundary layer known. 


Division the Boundary Layer into Regions 


The shear stress turbulent flow composed viscous and turbulent 
terms, and may written: 


(1) 


where the molecular viscosity and the eddy viscosity associated with 
the turbulent fluctuations. The usual approach the turbulent boundary layer 
terms the relative influence the molecular and eddy viscosities. 
Thus, very close the wall there the laminar sublayer which the molecu- 
lar viscosity predominates and which the velocity proportional the dis- 
tance from the wall. Further out the molecular viscosity has negligible ef- 
fect and the eddy viscosity governs the flow. between, there buffer zone 
which both viscosities must considered. Except the lowest Reynolds 
numbers, the laminar sublayer and the buffer zone are very small. One 
therefore left with single turbulent region and solutions the equations 
motion are sought which the molecular viscosity practically neglected. 

Usually single function used describe the entire turbulent velocity 
distribution. The well-known power law example such formula. 
Sometimes gives reasonable fit, but poor approximation. 
Superposition formulas developed Hudimoto 1935, and more recently 
Rotta(4) and Ross and Robertson also fail agree with the details the 
velocity profiles. Specifically, none these single functions admits the exist- 
ence inflections the profiles; and yet, such inflections are invariably 
found foi profiles close separation. 

Recognizing that the usual analyses not achieve sufficient agreement with 
the details measured velocity profiles, the physical factor was sought that 


Ross, Ph. Thesis, Ibid. 
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would account for the experimentally observed features turbulent boundary 
layers. This factor the relative proximity the wall. the new analysis, 
there superimposed the usual division according viscosity considera- 
tions division into three regions according the relative influence local 
wall conditions. These zones are approximately the inner ten per cent the 
boundary layer thickness, which the proximity the wall influences the 
flow; the outer seventy-five per cent, governed the spatial history the 
boundary layer; and, blending zone between the others which the flow 
influenced both the wall proximity and the spatial history. Fig. shows 
typical velocity profile and illustrates the superposition the two methods 
dividing the boundary layer. The inner region seen include the laminar 
sublayer and the buffer zone well inner turbulent part. 

understand the reasoning behind this method division essential 
realize the fundamental distinction between laminar and turbulent shear 
flows. laminar flows the property that gives rise the shear stress, name- 
the viscosity, exists throughout the fluid uniform manner. The laminar 
boundary layer exists because the physical boundary condition zero velocity 
solid surface differs from the potential theory solution for the main fluid 
flow. the case the turbulent boundary layer, the high intensity turbulence 
that supports the shear has its origin very close the wall and diffused into 
the fluid. Because the turbulence originates very close the wall, the re- 
gion very high velocity gradients, and diffuses out into the moving fluid, 
follows that the turbulence intensities and shears the outer part the bound- 
ary layer are governed the spatial history the flow while those near the 
surface are more function local conditions. This concept the space- 
history the flow was expressed Schultz-Grunow 1938 and again 
Dryden(7) 1946. was aware this important point early 
1935, when stated, "for each velocity profile develops some regular man- 
ner from the velocity profiles situated further upstream; i.e.....depends the 
previous history the portion the fluid This concept the 
diffusion turbulent stresses and the importance spatial history be- 
hind the division the boundary layer into the zones shown Fig. 

The properties the inner region are closely related the presence 
the wall. Here occurs the majority the conversion mean motion tur- 
bulence and turbulence heat. The production turbulence predominates 
over its dissipation. the outer region, both production and dissipation occur 
very much reduced rates, with dissipation predominating. The outer region 
also characterized "intermittency effect” that arises from the irregu- 
lar nature the outer boundary the turbulent region{9) 

Between the inner and outer regions there blending region that par- 
tially dependent local conditions and partially spatial history. This blend- 
ing region combines the features the other two regions and necessarily 
more complex than either the others. not confused with the "over- 
region. analyses equilibrium boundary layers Millikan 10) and 
von Mises(11) the velocity profile considered from the wall outward and 
from the outside in. also assumed that there overlap region close 
the wall which both functions are equally valid. the overlap region 
closer the wall than the blending region, the overlap condition can only 
used when the outer profile known through the blending region. Thus, the 
overlap concept has some apparent value when one dealing with equilibrium 
boundary layers, but value when dealing with the more general non- 
equilibrium case. 
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Velocity Functions 


the methods dimensional analysis are applied the problem turbu- 
lent-boundary-layer velocity profiles, using the concepts already developed, 
two The one for the inner profile: 


involves the wall-shear-stress the definition the friction velocity 


(3) 


the wall roughness the term y/k, and the local pressure gradient the 
last term. The relation for the outer profile: 


expresses the velocity deficiency terms the relative distance from the 
outer edge and the spatial history. (u, the free-stream velocity and 
the effective boundary layer thickness.) 

The inner function consists three parts: the laminar sublayer, the buffer 
zone, and the inner turbulent region. The first two these are usually quite 
thin. The inner turbulent region characterized logarithmic velocity 
profile, which can derived directly from dimensional analysis assuming 
that most important influence the eddy viscosity the distance from the 
wall. The equation, the form 


well ver measured velocity profiles. The Law the Wall proposed 
Ludwieg 13) 1949 states that for smooth walls the coefficients and 
the logarithmic inner profiles are practically universal constants indepen- 
dent pressure gradients, the Reynolds number, and the past history 
the flow. (The functional relation equation includes pressure gradient 
term, but this may shown small except for low Reynolds numbers 
and extreme pressure gradients.) The author has approximately 
one dozen experiments which both wall shear stress and velocity profiles 
were measured and has found appreciable difference the coefficients for 
the various types flows. For all practical purposes, the coefficients and 
are both 5.6, and the formula for the inner turbulent region is: 


For the outer turbulent region, simple formula which gives excellent 
fit most velocity profiles the three-halves power deficiency law: 


Ross, Opus Cit. 
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This formula, which valid for approximately the outer seventy-five per cent 
the boundary layer, similar that proposed Darcy, nearly century 
ago, for the central region pipe flow; hence the choice the letter for the 
coefficient. Its local value, which depends the spatial history the flow, 
varies from about 0.3 for flatplates high 1.3 near separation. 

The various assumptions that have been used derive the inner and outer 
velocity functions fail the blending region between these two regions. The 
physics this region not such that the velocity profile can described 
any simple functional relation; neither wall conditions nor spatial history can 
ignored. infinite number outer profiles may exist. For each outer 
profile variety inner profiles may exist, depending primarily the local 
Reynolds number. The velocity profile the blending region will depend upon 
the combination inner and outer profiles. For those cases for which the 
extrapolated inner and outer functions cross each other appreciable angle, 
the curve the blending region will include inflection point. Some the 
possible combinations are illustrated Fig. which shows clearly the in- 
flections that occur for high profiles. 

The division the boundary layer into inner and outer regions, superim- 
posed the normal division into viscous and turbulent zones, has resulted 
somewhat more complex picture the turbulent boundary layer than that 
usually considered. That this additional complexity worthwhile supported 
the simplifications that result the analysis each the regions. Fig. 
illustrates the division the boundary layer into regions and summarizes 
their properties. Along with the extents the regions and the types veloc- 
ity profiles, the general characteristics the shear stress profile and the 
Prandtl mixing length (14) are indicated. These are based experimental 
observations and are included because the application Prandtl's momen- 
tum-transfer theory later stage the analysis. Fig. shows 
set experimental profiles, measured Schubauer and The 
solid curves are the inner and outer profiles given equations The 
curves progress from equilibrium, flat-plate conditions practically separa- 
tion. The blending between the two portions clearly illustrated. 

The strongest argument for the present analysis that the simplest 
way treating the turbulent boundary layer consistent with known physical 
properties that leads useful results. The remainder the present paper 
devoted brief discussion the new method for calculating two-dimensional 
turbulent boundary layers that was developed from this analysis. 


Calculation Method 


When the turbulent boundary layer problem expressed its most general 
form there are four simultaneous equations, two which may non-linear 
differential equations and all which are interdependent. These equations are: 

the growth the momentum thickness function the wall shear 

stress, pressure gradient, and shape parameter; 

the shape parameter the velocity profile function the pressure 

distribution and the spatial history the flow; 

the wall shear stress function the shape the velocity profile 

and the momentum-thickness Reynolds number; and, 
the detailed velocity distribution function the shape parameter, 


604-5 


the boundary-layer thickness, and the local wall shear stress. 

Obviously, simultaneous solution these four equations impractical. What 
required approximate method whereby the four equations can solved 
separately and, necessary, second approximations made. The analysis pre- 
sented here makes such solution possible for wide variety two-dimen- 
sional boundary layers. 

The boundary layer any station calculated follows: 

The momentum thickness found function the pressure distribu- 
tion, through recently developed approximate integration the 
von momentum equation. 

The coefficient the three-halves power outer velocity distribution 
replaces the shape function this coefficient related 
the spatial history semi-empirical equation derived from 
Prandtl's momentum-transfer theory. This function then leads 
itself through empirical correlation. 

The wall-shear-stress coefficient related and the local momen- 
tum-thickness Reynolds number Ludwieg's "Law the Wall," and 
second empirical correlation. 
Finally, the velocity distribution the station calculated two parts. 
The outer part three-halves power deficiency curve having coeffi- 
cient The inner part logarithmic with its magnitude determined 
from the wall-shear-stress coefficient. smooth curve used blend 
the two regions. 
The method applies redistributing flows that proceed from separation 
towards flat-plate equilibrium; well the usual boundary layer ad- 
verse pressure gradient that goes from equilibrium towards separation. The 
occurrence separation also correlated with the coefficient 


Growth the Momentum Thickness 


The equation that almost invariably used for the growth the momentum 
thickness the von integral momentum equation: 


where the wall-shear-stress coefficient and the common shape 
parameter. equation, somewhat different form, was originally derived 
von 1921 for use with laminar boundary layers. Although 
extensively applied turbulent-boundary-layer problems this form, ac- 
tually for many turbulent flows the contributions the turbulent normal 
stresses must included additional term. The exact solution the cor- 
rected equation requires the simultaneous solution the equations for the wall- 
shear-stress coefficient and for the shape parameter. However, Ross and 
have recently developed approximate solution that may 
used directly. They distinguish between two cases: small pressure gra- 
dients, for which the boundary layer close the equilibrium condition, and 
large adverse pressure gradients. the first case, step-by-step solution 

equation utilized, with both and assumed constant for each step and 
each usually having the equilibrium, flat-plate value. For comparatively large 
pressure gradients, they develop simple algebraic, power expression: 
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(9) 


which gives good agreement with experimental results. The exponent, 
function the momentum-thickness Reynolds number the initial station, 
varying from about 5.0 transition (R, about 4.2 high momen- 
tum-thickness Reynolds number about shown Fig. Utilizing this 
result, the momentum thickness for any station may readily computed with- 
out first solving the other three equations, and without knowledge the shape 
parameter. Usually direct algebraic solution will possible for 
function only the potential velocity, pressure, distribution. The mo- 
mentum-thickness Reynolds number then given by: 


The wall shear stress will shown partially function this Reynolds 
number. 


Behavior the Shape Parameter 


most the published analysis methods single-parameter family 
velocity profiles assumed and differential equation written for the vari- 
ation the shape parameter with distance. This auxiliary equation must 
solved simultaneously with the momentum equation and the equation for the 
wall shear stress. the present development, the boundary layer divided 
into regions and single parameter characterizes the entire velocity profile. 
However, the inner profile, well the outer one, dependent the outer 
shape parameter that solution equation for the analogous 
step the calculation procedure. 

The fundamental physical property governing the dependence the 
outer velocity profile the spatial history the flow and its independence 
local wall conditions. The problem expressing this equation was at- 
tacked through use the momentum-transfer theory. Although this 
mixing-length theory has long been held contempt many theoreticians, 
recent 19) the statistical aspects turbulence have 
shown that good approximation for the calculation velocity profiles. 

Prandtl's theory gives differential equation for the velocity profile 
function the shear-stress and mixing-length profiles: 

typical shear stress shown Fig. The outer region essenti- 
ally linear function with intercepts effective (outer) wall shear 
stress, and the effective boundary-layer thickness, The mixing 
length, assumed constant the region outside the influence 
the wall. When the linear shear stress profile and constant mixing length are 
assumed, Prandtl's theory yields the 3/2-power velocity deficiency relation 
with 


(12) 
where the constant mixing length are defined Fig. 
simple algebraic manipulation one can create parameter that should 


function distance but not local conditions. This combination is: 
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ing flow having constant outer velocity, shown Fig. Here the data 


where u,, measured the initial station. can reasoned that the 
outer region turbulent boundary layer similar the shear flow 
wake jet, then vary with distance manner independent 
changes local wall conditions and pressure gradients. 

the side equation for the data published Schubauer and 
Klebanoff this experiment the flow had practically zero pressure gra- 
dient the 18-ft station and strong adverse pressure gradient after this 
point. (Separation occurred just ahead the 26-ft station.) The linear trend 
that was established the constant pressure region continued into the strong 
adverse pressure gradient. Another confirmation, this time for redistribut- 


Klebanoff and Diehl for the flow along smooth plate preceded two feet 
very rough surface shows the same linear trend. This confirms basic 
assumption the present development, that all turbulent flows obey the same 
general relations. Most the published analyses are particularly unsuccess- 
ful calculating the shape parameter for this type flow. the original 
study® several other examples and redistributing flows are pre- 
sented, and concluded that, for distance equal about fifty initial- 
boundary-layer thicknesses, the outer velocity profile governed rela- 
tion the form: 


provided the flow very closely two-dimensional. Following the linear region, 
the initial character the outer flow pretty much submerged the turbu- 
lence generated within this distance and new equilibrium established. 

The coefficient may calculated for those cases for which equili- 
brium boundary layer established region zero pressure gradient 
ahead region relatively strong pressure gradient, change wall 
conditions. has constant value, may calculated from the equilib- 
rium, constant-pressure, flow: 


For equilibrium profiles, all the quantities are related, and possible® 
relate simply the value the wall-shear-stress coefficient: 


(16) 


which turn function the momentum-thickness Reynolds number, 
and/or the relative roughness. 

cases where known, the value any station can calculated 
the thickness the boundary layer, known. The momentum equation 
yields that the ratio sufficient. the blending region were 
fully understood, might possible calculate this ratio; however, was 


Ross, Opus Cit. 
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necessary resort empirical correlation. 

the original study’, over one-hundred experiments were examined, and 
eighteen were selected for the empirical analysis. These experiments were 
selected the best representatives wide variety experimental condi- 
tions. being logical contention the present analysis that all types tur- 
bulent-boundary-layer flows must governed the same physical laws, 
least one experiment for each type two-dimensional configuration was cho- 
sen. The correlation with shown Fig. From this result, 
can calculated known. Actually, returning equation 14, itis 
that found. Equation may written the form: 


where the left-hand side the function given Fig, 10. Thus, from 
equation and Fig. 10, the value any station may found provided 
that the value and the initial conditions are The above method 
computing lays stress the initial spatial variation the outer boundary 
layer. follows that boundary layers with different initial spatial histories 
(different will behave quite differently when acted the same pres- 
sure distribution. 


The Wall Shear Stress 


equation the inner velocity profile given function the friction 
velocity. useful expression for this latter quantity results from writing the 
inner velocity equation for the momentum thickness distance. The equation 
for this fictitious effective velocity, simply: 


(18) 


can manipulated into the expression: 


This equation may solved numerically for function and 
The result found fit very closely the approximate expression: 


should noted that this formulation the wall friction relation for turbu- 
lent boundary follows directly from the universal Law the Wall pro- 
posed The friction analysis given here parallels that 
Ludwieg with the difference being that the momentum-thickness velocity will 
related the new outer parameter instead the usual shape param- 
eter The empirical correlation with shown Fig. The 
scatter the data somewhat greater than that for but consistent 
trend could found with any other variable and the scatter attributed less 


7.Ross, Opus Cit. 
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second order effects than experimental error. Fig. enables the rapid 
estimation the wall-shear-stress coefficient function the coefficient 
and the momentum-thickness Reynolds number, 


The Velocity Distribution 


The velocity distribution made the separate inner and outer profiles 
connected blending curve. Equation for the inner profile may written: 


Where all the required quantities have already been computed. Similarly, 
the outer profile given directly equation terms known quantities. 
smooth curve may used connect the curves blending region. 
Thus the present analysis leads readily the computation the velocity dis- 
tribution station, without requiring the simultaneous solution pair 
differential equations. 


Separation 


The method the present paper particularly well suited the prediction 
separation. The calculation procedure gives way determining and 
the occurrence separation may correlated with Fig. shows linear 
variation for the velocity the momentum thickness function This 
extrapolates zero equal about 1.45, establishing upper limit. The 
actual value will depend somewhat the Reynolds number and the down- 
stream flow conditions. good criterion for separation D=1.340.1. 


CONCLUSION 


the present paper, new approach turbulent-boundary-layer problems 
has been outlined. The analysis considers that most important factor govern- 
ing the velocity profile the relative proximity the wall. Consideration 
this factor has led the division the boundary layer into three regions. The 
properties the inner, wall, region depend local wall conditions and the 
outer region. The latter determined function the spatial history 
the flow. The resultant semi-empirical calculation method the only one de- 
veloped date that applies redistributing flows well the usual case 
flow adverse pressure gradient. The calculation procedure also 
mathematically more simple, the usual simultaneous differential equations 
are replaced algebraic equations. This means that analytic answers are 
now possible for many problems previously requiring numerical solution. The 
method also yields simple criterion for predicting separation. Its most ser- 
ious limitation that, the algebraic form, only applies for distance 
equal about fifty initial boundary layer thicknesses. 
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Fig. Typical Idealized Velocity Profiles. 
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Fig. Typical Shear-Stress Profile for Flow with Adverse Pressure 
Gradient. 
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Fig. 11. Correlation Effective Relative Velocity 
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c. Discussion of several papers, grouped by Divisions. 
4. Presented at the Atlanta (Ga.) Convention of the Society in February, 1954. 
e. Presented at the Atlantic City (N.J.) Convention in June, 1954. 


AMERICAN SOCIETY CIVIL ENGINEERS 


OFFICERS FOR 1955 


PRESIDENT 
WILLIAM ROY GLIDDEN 


VICE-PRESIDENTS 


Term expires October, 1955: Term expires October, 1956: 
ENOCH NEEDLES FRANK WEAVER 
MASON LOCKWOOD LOUIS HOWSON 


DIRECTORS 


Term expires October, 1955: Term expires October, 1956: Term expires October, 1957: 
CHARLES MOLINEAUX WILLIAM LaLONDE, JR. JEWELL GARRELTS 
MERCEL SHELTON OLIVER HARTWELL FREDERICK PAULSON 
LLOYD KNAPP SAMUEL MORRIS DON CORBETT 
GLENN HOLCOMB ERNEST CARLTON GRAHAM WILLOUGHBY 
FRANCIS DAWSON RAYMOND DAWSON LAWRENCE ELSENER 


PAST-PRESIDENTS 
Members the Board 


WALTER HUBER DANIEL TERRELL 


EXECUTIVE SECRETARY 
WILLIAM CAREY 


ASSISTANT SECRETARY TREASURER 
CHANDLER CHARLES TROUT 


ASSOCIATE SECRETARY ASSISTANT TREASURER 
WILLIAM WISELY CARLTON PROCTOR 


PROCEEDINGS THE SOCIETY 


HAROLD LARSEN 
Manager Technical Publications 


DEFOREST MATTESON, JR. PAUL PARISI 
Editor Technical Publications Assoc. Editor Technical Publications 


COMMITTEE PUBLICATIONS 
SAMUEL MORRIS, Chairman 
JEWELL GARRELTS, Vice-Chairman 
GLENN HOLCOMB OLIVER HARTWELL 
ERNEST CARLTON DON CORBETT 


